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Abstract: In the framework of the Standard Model we calculate the flavour non-universal 
correction for neutrino refraction in a neutrino background and verify a similar previous 
result for the case of ordinary-matter background. The dominant term arises at loop level 
and involves r leptons circulating in the loop. These ©(G^m 2 ) corrections to the tree-level 
potential provide the dominant refractive difference between and v T unless the medium 
contains /i or r leptons. Our results affect the flavour evolution of dense neutrino gases 
and may be of interest for collective three-flavour oscillations of supernova neutrinos. We 
spell out explicitly how these non-universal neutrino-neutrino interactions enter the flavour 
oscillation equations. 



Keywords: 



Neutrino Physics, Standard Model, Supernovae. 



Contents 



[ll Introduction |l| 

^ Calculation of the radiative corrections ||| 
[231 Tree level | 



2.2 Radiative correction for a background of charged fermions 



2.3 Extension to neutrinos as background fermions H 
|3L Including radiative corrections in the oscillation equations 11 
EL Conclusions 12 



1. Introduction 

The dispersion relation of any particle is modified by a background field or medium. For 
neutrinos this effect is extremely faint, yet the matter effect is crucial for neutrino oscilla- 
tions where the feeble refractive indices compete with the minuscule mass differences ]l|, |2| . 

Since neutral-current interactions among the three active neutrinos are flavour univer- 
sal, they can be usually ignored for neutrino oscillation phenomenology, where only the 
relative energy shift between flavours is important. Moreover, energy scales of interest are 
well below the weak boson mass Myy, so for low-energy neutrinos in ordinary matter (no 
/x or t leptons) the potential-energy shift between flavours can be calculated simply from 
the effective charged-current Hamiltonian density 

T-Cmsw = -7=[e7a(l - 7sKPe7 Q (1 - 75) e] = -^[v e y a (l - l5)^e][ej a (l - 75 )e] , (1.1) 

where the second equality follows from the usual Fierz rearrangement. One thus finds the 
well-known tree- level energy splitting for electron neutrinos with respect to or v T , 

AE UeU ^ T = E Ve - E^ T = V2G F n e , (1.2) 

where n e is understood as net density of electrons minus positrons. In lowest order of 
perturbation theory, the f M and v T dispersion relations coincide and AE Ufil/T = 0. 

Equation ( |1.2| ) provides a sufficient description of the matter effect in most practical 
cases, most notably solar neutrinos, but care must be taken in some circumstances where 
either the leading term vanishes (n e ~ 0) or subtle three-flavour effects come into play. 
For example, it has long been appreciated that radiative corrections break the v^-Vt de- 
generacy in an ordinary medium because the charged- lepton mass (m^ or m T ) in the loops 
distinguishes between flavours [||. In the low-energy limit, the loop-induced v^-Vj- energy 
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shift has the same structure as Eq. ( |1.2| ) and thus can be represented by an effective net 
tau-lepton density in the form 

AE VTU)t = E Ut -E^ = V2 G f nf . (1.3) 

Botella, Lim and Marciano || studied explicitly the case of an ordinary electrically neutral 
medium and found the equivalent of 



cff 3 



n B , (1.4) 



where ng = n p + n n is the baryon density and Y n = n n /riB the neutron/baryon fraction, 
the remainder being protons. We have introduced the small parameter 

c = ^ = 264 x 10 <L5) 

that will appear in all of our results. Here and in the following we neglect m e , rria, the 
first-generation quark masses and the neutrino energy relative to m T and Mw- 

The quantity Y n /3 in Eq. ( |1.4| ) is at most a 5% correction. Neglecting it provides 

nf ~ 2.6 x 10~ 5 n B . (1.6) 

The only known example where this small density can be of practical interest is the case of 
supernova (SN) neutrinos propagating in a dense environment after leaving the neutrino 
sphere of the collapsed core. At large distances they encounter two resonances where 
the usual matter effect causes a level crossing with the "atmospheric" and "solar" mass 
differences, respectively Q. Much closer to the core the radiative effect will cause similar 
resonances ||. However, the u^, u^, v T and v T fluxes and spectra are probably equal, so 
one may think that v^—Vr transformations play no role. 

This picture has radically changed with the insight that typically flavour oscillations 
do occur close to the neutrino sphere. The refractive effect of neutrinos on neutrinos is not 
fully captured by a refractive energy shift because if neutrinos oscillate, the same applies 
to the neutrino background, providing a "flavour off-diagonal refractive index" || f?j]. In 
other words, the neutrino oscillation equations become intrinsically nonlinear, leading to 
collective forms of oscillation |||, ||], a subject of intense recent investigation (see the recent 
review Jul and references therein) . Since collective oscillations operate close to the neutrino 



sphere, here the v^,~v T resonance can become important in a three-flavour treatment [11 



12]. On the other hand, collective oscillations are suppressed if the matter effect is too large, 



typically when ng > n u [13]. The SN matter profile and neutrino flux vary with time and 
depend strongly on the properties of the progenitor star, so depending on circumstances n v 
can be larger or smaller than ub in the collective oscillation region a few hundred kilometers 
above the neutrino sphere. When n v > jib and non-linear effects become important, it is 
reasonable to expect that the leading effect breaking the degeneracy between i/„ and v T 
is actually provided by loop-induced effects due to the neutrino background itself. In this 
article we fill a gap in the literature and calculate this non-universal radiative correction 
for neutrino refraction in a neutrino background. 
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Our calculation is described in Sec. |2| where we begin with a background of charged 
fermions, confirming the original result shown in Eq. ( |1.4|) . From this starting point we in- 
clude neutrinos as a background and derive our main result, an effective neutrino-neutrino 
Hamiltonian that includes m T effects in the loop. In Sec. |3| we spell out how to include our 
result in the non-linear flavour oscillation equations before concluding in Sec. |j. 



2. Calculation of the radiative corrections 

The energy shift between different neutrino species in a given medium is proportional to 



the difference of forward scattering amplitudes [14, 15 1. Since we deal with energies well 
below the electroweak scale, the results of the one-loop calculation within the electroweak 
theory can be expressed in terms of an effective four-fermion Hamiltonian. Generically, the 
quantum effects manifest themselves in the low-energy theory as corrections to the coeffi- 
cients of the operators. (In addition they could induce higher-order operators mediating 
processes forbidden in the original theory.) The corrections, however, need not respect the 
original symmetries of the interaction. The case at hand makes no exception: we shall 
see by explicit calculation that at one loop, although the Lorentz structure of the tree- 
level current-current interaction is preserved, flavour universality of the "neutral currents" 
breaks down. 

2.1 Tree level 

To set up our problem and fix the notation, we recall that neutrinos interact with matter 
via the charged and neutral weak-current Lagrangian densities 

C CC = --^=(J% c W +a + J^W~ a ) , (2.1) 



2\/2 

9 
2c w 



£nc = -t^-^ C ^- (2.2) 



Here, g denotes the SU(2) weak coupling, c w = y^l — s\ is the cosine of the weak mixing 
angle, and and Z represent the charged and neutral gauge bosons with masses Myy 
and Mz = M\y/c w , respectively. These parameters are related to the Fermi constant Gp 
and the fine structure constant a by 

_ V2g 2 _ to 

F "^"w (2 - 3) 

In terms of the weak-interaction eigenstates, the charged current is 

J a° = Yl ^Ta(l " 75^ + Yl " TfeK (2-4) 
I i 

while the neutral current is 

J « C = Yl [^{CV ~ CAl5>l + ha{c l V ~ C l Alh )l\ + 
I 

+ iZ ["^a(cv - CAl5)ui + di7a(cy - c d A ^)di . (2.5) 



-3- 



The indices I and i run over the lepton and quark generations, respectively. The vector 
and axial couplings of a fermion / = u, e, u or d are 

4 = T 3/ -2Q / 4, c{ = T 3f , (2.6) 

where Qf = {0, — 1,|,— |} is the electric charge and T%f = \, —\] is the third 

component of the weak isospin. In particular c v = c A = 1/2 for neutrinos. We need to 
calculate the shift in the potential of a neutrino v\ = u^, v T interacting with a medium 
consisting of ordinary matter, i.e. fermions / = u,d, e. This interaction is described by the 
effective neutral-current Hamiltonian density 

H nf = ^§fo7«(l - lMh a {c v - c Al5 )f] . (2.7) 

Charged-current interactions can be ignored since we consider only the potential-energy 
difference between [i- and r-neutrinos, which is a purely radiative effect, and the medium 
is free from [i- and r-leptons. 

Assuming homogeneity, the potential in terms of the Dirac spinors of the (left-handed) 
neutrino u u {p) and of the fermion Uf(k, s) is 

V «f = %^j d3kd3xi? H k ' s ) MP)7«(1 " J 5 )uv( P )][uf(k,sh a (c v - c A ^)u f {Ks)] 

(2.8) 

where the spatial integral is over the normalization volume (containing one neutrino), 
Ff(k,s) is the momentum distribution of the fermion background (normalized to 1), and 
s denotes the fermion polarization. The matrix elements are calculated between identical 
initial and final states, both in momentum and spin, i.e. they correspond to the forward 
scattering amplitude. This leads to [16| 

v -t = S E/ d 3 kf>(k,s) (§?) [s/(M)7«(4 -<&*)«/(*.< 

- V2 F E / *kF,M (|) (n /v k °- E j™^ , ( 2 . 9) 
where is the background fermion mass and s a its polarization four vector, such that 

8 ■ k = 0, S 2 = -1. 

When the medium consists of unpolarized ordinary-matter fermions (/ = e,u,d), the 
term linear in s a averages to zero. Thus the neutrino potential is independent of the axial 
coupling c A and reads 

V Ulf = V2G F n f c f v I d 3 kF / (k, S )(l-v p -v k )=V2G F n/4, (2.10) 

s 

where v p and Vk are the neutrino and fermion velocities, with momentum p and k, re- 
spectively. In the last equality we also assumed isotropy of the background distribution. 

As far as the tree- level neutral-current interaction in Eq. (2.7) is concerned, Vn, v T , 
and v e receive a common potential shift described by the flavour-independent result of 
Eq. ( 2.10 ), which is irrelevant for neutrino oscillations. Note also that, in a neutral medium, 
the sum over all species cancels the terms proportional to the electric charge Qf in Eq. fl2.G| ) . 
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Figure 1: At tree level, neutrino forward scattering in a charged-fermion background (/ = u, d, e) 
proceeds through t-channcl exchange of Z bosons. 



2.2 Radiative correction for a background of charged fermions 

The effect we are interested in is the flavour-dependent part of the loop contribution to 
the neutral-current couplings Cy, A in Eq. (^7j). Technically, the matching of the effective 
Hamiltonian coefficients will be obtained by computing the v\-f forward scattering ampli- 
tude. Note also that the momenta of the scattering particles will be set to zero, since their 
typical values in SNe are very small compared to the relevant mass scales m T and M\\r. In 
this approximation all one-loop integrals reduce to tadpoles. 

At tree level, v\ forward scattering in a charged-fermion background is described by 
the Feynman diagram in Fig. The corresponding amplitude is 

• 2 ■ 2 

A " lf = 5J^ 3 ^*7m(4 - ^75)^2] = ^rW3ic£ , (2.11) 

where ojl = (1 — 7s)/2 is the left-handed projector and we have introduced the currents 

C{j = ua^{c f v - c f A -f 5 )uj , Wij = Cfj = Ui ^uj l uj , (2.12) 

where the Lorentz indices and their contraction are implicitly understood. At one loop, 
neutrino-flavour dependent contributions result only from diagrams where a charged gauge 
boson (W^) or would-be-Goldstone boson ((p^) couples to a neutrino line, giving rise to 
a virtual lepton. In practice we will restrict ourselves to this type of diagrams and we 
will retain only terms that depend on the mass of the lepton and are thus sensitive to 
its flavour. Flavour-independent contributions will be systematically neglected. Note that 
here and throughout we work in the 't Hooft-Feynman gauge, in which = Myy- 

As we will see, it is convenient to classify the loop diagrams according to their topology 
and to treat box and non-box diagrams separately. The latter represent the corrections to 
the viV[Z vertex and comprise the vertex diagrams of Figs. ||-|3] as well as the selfenergy 
insertions of Fig. ||. Note that all diagrams of non-box type are independent of the nature 
of the background fermions / in the sense that the tree-level ffZ vertex factorizes. Box 
diagrams can involve two charged bosons (WW , W<f> or cjxp) and are suppressed unless 
both of them are W bosons. Depending on the isospin nature of the background fermions, 
different box topologies contribute (Fig. ||): in presence of a down- type background (/ = 
e, d) only ladder boxes appear, while crossed boxes contribute only if the background 
fermions are of up type (/ = u and, for later application, also / = u). 
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Figure 2: Flavour-dependent vertex corrections involving Zll interactions. 




Figure 3: Flavour-dependent vertex corrections involving ZW + W~ , Z4> + 4>~ , and ZW + <fi~ inter- 
actions. In the latter case (rightmost diagram), an additional contribution with W <-> 4> must be 
taken into account. 




Figure 4: Flavour-dependent selfenergy insertions in the incoming neutrino line. Analogous di- 
agrams for the outgoing neutrino must be considered. These diagrams have to be understood as 
on-shell renormalization factors for the neutrino wave functions. 



Apart from flavour-independent additive terms, the relevant amplitudes for the dia- 
grams in Fig. H are 



+ 



m. 



4(4tt) 2 M^ 



2 + * 



-flog 



raj 



mic f 



42 i 



mr 



4(4vr)2 m; 



1 i A* 
A+ 2 +1 ° g M? 



2 1 



W 



(2.13) 
(2.14) 
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Figure 5: Flavour-dependent box corrections. The ladder box diagrams (left) contribute only for 
down-type (/ = e, d) background fermions, whereas the crossed box diagrams (right) contribute 
only for up-type (/ = u, v) fermions. Box diagrams involving ^ exchange are suppressed. 



For the diagrams in Fig. |3| one finds 
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4(4vr) 


2 ^ 



A 1 ^ 
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(2.15) 
(2.16) 
(2.17) 



where in the last equation we included a factor two corresponding to a second diagram 
with W <-» 0. Finally, the relevant amplitude for the diagrams in Fig. is 



A' 



m, 



4(4vr) 2 M 4 V8 4 



1, M 
i log- 



W3lC{ 2 



(2.18) 



where, again, we included a factor two for the correction to the second neutrino. The 
expression A = — 7_e + log 4-7T represents ultraviolet poles within dimensional regular- 
ization and /i is the corresponding mass scale. The complete non-box contribution, i.e. the 
sum of the amplitudes Eqs. (2.13)~( ^.18 ), is free from ultraviolet divergences and reads 



234 4(4^) 2 M 4 



w 



3 + log 



m 



2 1 



M w 



W31C 



42 



(2.19) 



For the box diagrams in Fig. 5, neglecting contributions of order mf/Mw, one obtains 



A VlU 



A 



m 



' " (4tt) 2 M 4 



1 ig" 



in 7 



4(4vr) 2 M 4 , 



w 
log 



log 



m l 



+ 1 



W31W42 , 



+ 1 



W31W42 



(2.20) 
(2.21) 



The results for down- and up-type fermions correspond to ladder and crossed box topolo- 
gies, respectively. 

These corrections can be described as flavour-dependent shifts, 



J 

-V/A 



J 



r -> _|_ A r v 



(2.22) 
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of fermion couplings in the effective Hamiltonian density 7i Ul f of Eq. (2.7). The shifts 
induced in the vector and axial couplings of down- and up-type fermions are 



Ac 



Ac 



vid _ 


9 2 


2 

mf 


V/A ~ 


(4vr) 2 


M w 






2 


UlU _ 


9 2 


mf 


V/A ~ 


(4vr) 2 





3d 



V/A 



+ 2 + 



4 M + 2)log^ 



2 1 



-V/A 



W 



(2.23) 
(2.24) 



Here the terms proportional to Cy^ A originate from non-box diagrams and the remnant is 
due to boxes. The relevant shift for the electron component of the background is obtained 

flp, TL U = 2 Tip + Tl n , fid = 



by replacing d — > e in Eq. (|2.23[) . 

Taking into account that in neutral, ordinary matter n e 



n p + 2 n n and using n b 



n p + n r 



V2G F 7" A C y rf n f 



one obtains the Vn~ v T energy shift 



/ 



9' 



mz 



log 



ml 



1 + 



n B , (2.25) 



equivalent to Eq. (|0|) and in full agreement with Ref. ||. 
2.3 Extension to neutrinos as background fermions 

Now we turn to our main case of interest, namely when the propagating neutrinos {y{) 
interact with a neutrino background (/ = vy\ At tree level the relevant interactions are 
described by the four-neutrino neutral-current effective Hamiltonian 

G F \. — v 



(2.26) 



which is of the same form as Eq. ( |2.7| ) upon identification of the fermion couplings with 
Cy = c A 1 ' =1/2 and noting that each term with I ^ V enters twice in the sum. Note also 



that Eq. fl2.26j ) has a U{3) symmetry in flavour space. 

The general form of the one-loop effective Hamiltonian is 

H vv = Y, U ^v = % zD 1 + kUIUi ') WlcVM] fo'-fw] • (2-27) 
1,1' ^2, l v 

Again only flavour-dependent contributions k" 1 " 1 ' are taken into account. In practice we 
consider only terms involving mj/M w and/or mf,/M w . Moreover, the electron and the 
muon masses are neglected, i.e. we use 



m i 



5h 



m T 



This implies 



0, 



(2.28) 



(2.29) 



and one-loop terms proportional to m^/M w contribute only if either v\ or vy is a v T . 

We first consider the correction n UrUfi = k v $ Vt with /3 / t. This case can be easily 
related to the charged-fermion background results of Sec. |2.2| . To this end we can regard 
v T as a neutrino that propagates in a background with / = vq. The one- loop diagrams 



that give rise to terms of order m^/M^y are exactly the same as in Sec. 2.2 
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In principle, in addition to the vertices and selfenergies of Figs, |2j-||, which contribute 
to the v T v T Z vertex, diagrams corresponding to the vpvpZ vertex should be included. 
However, for (3 ^ r, such diagrams are free from m^/MrL contributions and can thus 
be neglected. 

Concerning the box diagrams of Fig. ||, the terms of order mj^/M^ that we have 
neglected in Sec. [T^ correspond to terms of order rnl/Myy and remain negligible. 



This implies that k VtV i 3 can be obtained from Eq. (2.24) for up- type fermions by simply 
replacing u — ► vp. By comparing the corresponding Hamiltonians Eqs. ( p. 7] ) and (2.27) we 
easily see that the relevant relation is 

Ac UtUi3 

K VrV g= _m =2Ac ^_ (2 _ 30) 

C V/A 

Inserting the explicit neutrino couplings we obtain 

K v T vp — k V0Vt — _o 9 m r _ _ ( Q Ol \ 



where e was defined in Eq. Q1.5|). 




Figure 6: At tree level, neutrino forward scattering in a neutrino background proceeds through 
i-channel exchange of Z bosons (left diagram), plus it-channel exchange in the case of identical 
flavours only (right diagram). 

Finally we consider the equal-flavour corrections k VtVt where the following differences 
must be taken into account: 



All scattering amplitudes are a factor two larger because in addition to the t-channel 
Z-exchange diagram that contributes to the v\ ^ vy interactions (left diagram in 
Fig. |6|), also a u channel opens (right diagram in Fig. |6|). Moreover, due to a Fierz 
identity, the amplitudes of the t- and u-channel diagrams are identical. In particular, 
the tree-level v T 



amplitude reads [cf. Eq. (|2.11[) 1 



A? Ut = 2 x A\ 



2 



W31W42 • 



(2.32) 
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In addition to the vertices and selfenergies of Figs. |2|-||, similar diagrams must be 
included, that describe the correction to the ffZ vertex with / = v T . In practice 
these two sets of "non-box" diagrams provide exactly the same correction. Includ- 
ing a further factor of two due to the presence of u-channel contributions we have 
[cf. Eq. ([HQ}] 



AV T V T 

• /i 234 



2 X 2 X 



A"Tf 

• /i 234 



1 9 



mr 



(4vr 



3 + log 



m i 



w 31 w 42 



(2.33) 



Finally the contribution of the box diagrams of Fig. |^ cannot be inferred from the cal- 
culation performed for the matter-background case since the terms of order m 2 ,/M' 2 y 
cannot be neglected when / = vy = v T and /' = V = r. In this case, including 
again a factor two due to u-channel contributions, the ladder box diagram yields 



[cf. Eq. (Ht)] 



AV t Vt 

•^5 



*5 



mz 



2 (4vr)2 



3 + 21og 



W 31 W 42 



(2.34) 



Combining the above results we obtain 



g 2 m 2 . 



(4tt) 2 M, 



w 



-2 e - 



(2.35) 



where e was defined in Eq. ( |1.5| ). 

In the context of neutrino oscillations (see Sec. |3|) it proves useful to rewrite our results 
in a more compact notation: 

H VV = ^L{ W - eT) 7 W] W - eJ) lx oj L v] + | [i>T^uj L v} [vT^lv] } , (2.36) 



V2 

where v is now a vector in flavor space and in the weak-interaction basis 

T = diag(0, 0, 1) , 



(2.37) 



i.e. the matrix T is a projector on the r direction in flavour space. Equation ( [2.36 ) holds 
at leading order in e. Note that the renormalization of the v T -v T coupling implied by the 
first term in Eq. ( |2.36| ) is 1 — 2e and has to be summed to the second term to provide 
the correct result 1 — 3e/2. In other words, not only is the original U(3) symmetry of the 
tree-level Hamiltonian Eq. ( [2.26| ) broken, but a simple description by renormalizing the 
coefficients Cy^ A is not possible. 

For comparison with the refractive energy shifts caused by a background of charged 
fermions, we recall that at tree level a neutrino v\ in a homogeneous and isotropic bath of 
v\i experiences 

(2.38) 



V Vl v v =V2G F {l + 8 w )n Vl ,. 



Of course, the usual ultra-relativistic approximation was used where neutrinos are purely 
left-handed and n V l is once more a net density of neutrinos minus antineutrinos. If neither 
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I nor I 1 are the r flavour, radiative corrections are negligible at our level of approximation. 
For I = r and I' = (3 7^ r our results imply 

V VrVp = V2G F (1 - e) n U0 (2.39) 

and analogously under the exchange of r *-* (3. Finally, for I = V = r we find 

V VtUt = V2G F (2 - 3e) n„ T . (2.40) 

However, in the context of neutrino oscillations these energy shifts do not provide a useful 
understanding of the role of neutrino-neutrino refraction because the background neutrinos 
themselves are, in general, coherent superpositions of weak interaction states. 

3. Including radiative corrections in the oscillation equations 

To see how neutrino-neutrino refraction influences flavour oscillations it is most econom- 
ical to represent the ensemble by flavour-space matrices of occupation numbers Qij(p) = 
(a\(p)aj(p)) where i and j are flavour indices, a} and a are creation and annihilation op- 
erators and (• • •) is an expectation value for the ensemble |l7|]. The diagonal elements of 
these flavour-space matrices are the usual occupation numbers, whereas the off-diagonal el- 
ements encode relative phases between the different flavour components. The simultaneous 
treatment of neutrinos and antineutrinos becomes particularly transparent if one defines 
Qijiv) = (^j(p)oi(p)) with "crossed" flavour indices where all over-barred quantities refer 
to antiparticles The equations of motion (EOM) for the flavour evolution are 

W P = [%,Qp] and i£ p = [Q p , g p ] , (3.1) 

where [•,•] is a commutator and the "matrices of oscillation frequencies" ft p and Q p play 
the role of Hamiltonian operators in flavour space for the evolution of mode p. 

In an isotropic ensemble, vacuum masses, matter, and neutrinos provide at tree level 
the usual terms |7|] 

% = +^ + V2G F (Nj + ISy , Q p = + V2G F (N z + N„) , (3.2) 

where M 2 is the neutrino mass-squared matrix and E = |p| for ultrarelativistic neutrinos. 
In the weak interaction basis, M 2 = UDU^ where D = diag(mf , m|, m|) and U is the 
usual mixing matrix. The matrix of number densities for charged leptons is defined 
in the weak-interaction basis with diagonal elements ni = rii- — rii+ . The corresponding 
neutrino matrix is 

If the medium is not isotropic, as in the case for neutrinos streaming from a SN core, the 
usual factor (1 — v q • v p ) must be included under the integral and the matrix then also 
depends on p. This structure is identical in all that follows, so without loss of generality 
we can study the radiative modification assuming an isotropic medium. 
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To see how our radiative corrections modify the structure of N^ it is easiest to begin 
with a hypothetical case of non-universal neutrino interactions of the form 

Hq = -j=(uGj uj l u)(uGj\uj l u) , (3.4) 

where G is a matrix in flavour space of dimensionless coupling constants. The tree- level 
standard-model case is G = 1 (unit matrix in flavour space). It was previously shown that 
the oscillation equation is to be modified in the form @ 

N„ - Nf (G) = j-^{G(ec l - g q ) G + GTr[(f? q - £ q )G]} . (3.5) 

For G = 1 this expression is identical to N u up to an additive term that is proportional to 
the unit matrix and to the net total neutrino density 



*' q Tr(£ q - £ q ) . (3.6) 



(2tt) 3 

This quantity is conserved under oscillations, so and Nf (1) play identical roles and 
actually we could subtract Eq. (3.6) from Eq. ( |3.5[ ) to achieve H v = Nf (1). 

Our calculation of the radiative correction of neutrino-neutrino interactions reveals 
that the flavour-sensitive difference of the scattering amplitudes can be expressed, at leading 



order and in the "large m T limit," in the form of Eq. ( p.36| ). Therefore, radiative corrections 
effectively enter with two different matrices of non-universal couplings. In the oscillation 
equation we thus need to substitute 

N^Nf(l-eT) + |Nf(T), (3.7) 

where the projection matrix T was defined in Eq. (2.37). In the weak- interaction basis this 
is explicitly to leading order in e 

U e J 

N^N^-el () N£ T | . (3.8) 

K e + N ^ + 2N7 

where we have neglected a term proportional to the unit matrix that is irrelevant for 
oscillations. 




4. Conclusions 

We have calculated the flavour non-universal correction for neutrino-neutrino refraction in 
the framework of the Standard Model. The dominant term arises at loop level and involves 
r leptons circulating in the loop. In the course of our derivation we have reproduced 
a similar term for a background medium of charged fermions that had been calculated 
previously by Botella, Lim and Marciano ||. One novel feature of our result is that the 
radiatively corrected effective neutrino-neutrino Hamiltonian can not be expressed in terms 
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of renormalized tree-level coupling constants. This different structure derives from the box 
diagrams with the exchange of two W bosons. 

The impact of neutrino-neutrino refraction on collective neutrino oscillations is not eas- 
ily assessed by comparing refractive energy shifts relative to those caused by an ordinary- 
matter background. The recent torrent of activities concerning collective SN neutrino 
oscillations was essentially triggered by the insight that the nonlinear nature of the equa- 
tions allows for large collective effects even if the ordinary-matter background causes larger 
energy shifts than the neutrino background. This is traced back to the phenomenon that 
the refractive-index matrix caused by ordinary matter and by background neutrinos are 
usually not diagonal in the same basis because of neutrino flavor oscillations. 

Our radiative corrections provide non-universal neutrino-neutrino interactions similar 
to, but not identical with, the case of non-standard interactions studied in Ref. 18], sug- 
gesting the possibility of interesting collective oscillation phenomena. However, a realistic 
assessment requires a dedicated study including the impact of a dense ordinary-matter 
background and multi-angle effects. 

We also note that in certain supersymmetric scenarios the neutrino-fermion radiative 
corrections could be enhanced up to two order of magnitudes with respect to the standard 
case [IS, 20]. Similar enhancements might be expected for neutrino-neutrino refraction. 
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